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ANOMALOUS  COMPONENTS  OF  SUPERCOOLED  HATEK  EXPANSIVITY,  COMPRESSIBILITY 

AND  HEAT  CAPACITY  (C  and  C  )  PEON  BINARY  FOBMAMIDE  ♦  HATEK  SOLUTION  8TUDIES 
P  ▼ 


M.  Oguni  and  C.  A.  Angel 1 
Department  of  Chemistry 
Purdue  University 
Heat  Lafayette,  Indiana  47907 


Abstract 

Following  heat  capacity  studies  of  Njl^  ♦  HjO  and  H202  ♦  HjO  solutions  fron 
which  an  anomalous  component  of  the  pure  water  behavior  could  be  extracted  by 
extrapolation,  we  extend  this  type  of  study  to  a  system  which  has  the  chemical 
stability  needed  to  permit  expansivity  and  compressibility  measurements  as 
well.  Data  accurate  to  ♦  2%  on  each  of  these  properties  aa  well  as  of  the  heat 
capacity  are  reported.  The  expansivity  data  support  almost  quantitatively  an 
earlier  speculative  separation  of  the  bulk  and  supercooled  water  expansivity 

into  a  "normal"  part  and  an  "anomalous"  pert,  the  latter  part  fitting  a 

.  -y 

critical  law  amnou  “  A(T/TB  -  1)  with  exponent  T  •  1.0.  According  to  the 
present  analysis,  the  excess  expansivity  which  is  always  negstive,  yields  I(  in 
the  rsnge  225-228,  y  in  the  rsnge  1.28-1.0,  depending  on  the  choice  of 
bsckground  extrspolstion  function.  The  bsckground  contribution  to  the  best 
cspscity  obteined  from  the  present  work  is  intermediste  in  chsrscter  to  thst 
from  the  previous  two  systems  sad  lends  to  similsr  equstion  psrsmeters.  The 
bsckground  contribution  to  the  compressibility  on  the  other  bend  is  very 
different  from  that  speculated  earlier  by  Kanno  and  Angel 1  and  approximately 
verified  by  recent  studies  on  ethanol-water  solutions  by  Coade  et  al/^  The 


background  component  from  the  present  analysis  ia  ^  50X  larger,  with  the  result 
that  the  excess  or  anomalous  component,  at  least  when  values  shove  0°C  are 
included  in  the  analysis,  cannot  be  sensibly  fitted  to  the  critical  point 
equation.  The  possible  origin  and  interest  content  of  these  differences  is 
discussed.  Combination  of  the  new  thermodynamic  data  permit  estimation  of  Cy 
values  for  the  solution,  and  by  extrapolation,  a  "normal"  Cy  component  for 
water.  The  "anomalous"  component  of  Cy  for  pure  water  obtained  by  difference 
has  the  form  of  a  Shottky  anomaly,  in  contrast  with  the  corresponding  Cp 
component  which  diverges. 


ANOMALOUS  COMPONENTS  OP  SUPERCOOLED  WATER  EXPANSIVITY,  COMPRESSIBILITY 
AND  HEAT  CAPACITY  FROM  BINARY  FORMAMIDE  ♦  WATER  SOLUTION  STUDIES 


M.  Oguni  and  C.  A.  Angell 
Department  of  Chemistry 
Purdue  University 
West  Lafayette,  Indiana  47907 

Introduction 

Recently,  in  an  attempt  to  elucidate  the  origin  of  the  many  anomalies  found 

in  the  behavior  of  supercooled  water,  we  reported  investigations  of  the  heat 

capacity  of  binary  solutions  of  water  with  the  water-like  second  components, 

(2) 

hydras ine  and  hydrogen  peroxide.  The  object  of  the  studies  was  to  obtain 
data  from  binary  aolution  regions  in  which  the  anomalies  characteristic  of  pure 
water  had  been  eliminated,  so  that  extrapolations  back  to  pure  water  would 
yield  an  estimate  for  the  "normal"  component  of  the  total  water  heat  capacity, 
thereby  permitting  the  magnitude  of  the  "anomalous"  component  to  be  extracted. 

Although  the  "normal"  component  obtained  from  the  two  binary  systems  was 
not  unique,  the  "anomalous"  component  thereby  obtained  showed  a  temperature 
dependence  which  conformed  well  to  a  critical  point  type  of  equation 

Cj  -  A(T/Tf  -  1)"Y  (1? 

with  a  singularity  temperature,  T^,  of  225  ♦  1.5  K  and  "critical"  exponents  Y 
either  aide  ef  the  value  of  unity,  depending  en  the  background  component  used. 

The  singular  temperature  is  in  good  agreement  (approximately  5°  below)  the 

1 

i 

1 

j 
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value  assigned  to  Tg  from  the  fits  of  total  supercooled  water  properties  to  the 

sane  equation  in  previous  analyses  of  the  low  temperature  anomalies. 

Significantly,  it  lies  just  a  few  degrees  below  the  lowest  temperature  to  which 

water  in  the  purest  condition  and  smallest  particles  yet  studied  can  be  cooled 

without  crystallisation  occurring,  a  finding  for  which  different 

(3-5) 

interpretations  have  been  advanced. 

It  is  clearly  desirable  to  determine  the  extent  to  which  these  observations 
can  be  confirmed  by  analysis  of  similarly  separated  anomalous  components  for 
the  other  intensive  thermodynamic  variables,  expansivity  and  compressibility. 
Unfortunately  neither  of  the  previously  studied  binary  systems  is  easily 
investigated  in  the  supercooled  regime  using  simple  dilatometric  technique 
developed  earlier  for  supercooled  water  studies,  because  the  second  component 
is  inherently  unstable  and  the  solutions  decompose  in  the  presence  of  the 
mercury  marker  bead  uaed  in  the  volumetric  analysis.  We  have,  therefore, 
sought  a  third  system  in  which  this  problem  would  not  arise  and  in  which  all 
three  thermodynamic  properties  could  be  studied  for  mutual  comparison. 

Hydroxyl amine,  HH^OH,  is  an  appealing  possibility  because  it  has  the  same 
lone  pair:proton  ratio  as  water,  in  constrast  with  (4:2)  and  H^H^  (2:4). 
However  HH^OH,  like  H^  and  ,  lacks  stability  and  therefore  was  not 
chosen.  After  consideration  of  several  possibilities  we  decided  on  formamide 
HHgCHO  which,  in  addition  to  being  stable,  seems  to  form  rather  ideal  solutions 
with  water.  The  m.p.  -5°C,  b.p.  125°C  and  dielectric  constant  c  "  55,  are  all 
comparable  with  water.  Ho  compounds  form  in  this  system  unlike  the  case  with 
H.O.  and  H„H. .  Both  ends  of  the  H1.CH0  molecule  participate  in  H-bonding  and 
the  lone  pair:proton  ratio  3:2  is  close  to  that  of  water.  Their  distribution, 
also,  is  not  completely  incompatible  with  the  open  bonding  of  water  since  the 
partial  molar  volume  of  water  in  the  binary  system  is  only  slightly  less  than 
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that  of  pure  HjO  across  the  entire  binary  system.^ 

He  have  succeeded  in  obtaining  data  of  good  quality  for  each  of  the  three 
thermodynamic  properties  for  solutions  at  20  mol  Z  intervals  across  the 
composition  variable.  This  paper  reports  the  results  of  these  aeasureaents, 
and  their  analysis.  He  have  also  studied  a  transport  property,  the  viscosity, 
but  will  report  the  results  of  that  work  separately.  Our  results  for  the 
coapressibility  aake  interesting  coaparison  with  the  current  findings  of  Conde 
et  al  for  ethanol  +  water  solutions  obtained  by  light  scattering 
aeasureaents. 

Experimental 

The  HH^CHO  used  was  Reagent  grade.  DSC  fusion  analysis  for  water  content 
(described  in  our  previous  paper)  showed  the  presence  of  0.25  mol  Z  H^O,  which 
was  allowed  for  in  all  solution  preparations.  BgO  used  was  deionized  triply 
distilled,  and  solutions  were  prepared  by  simple  weighing  procedures. 

The  heat  capacity  of  the  solutions  was  measured  over  the  temperature  range 
froa  25°C  down  to  a  low  teaperature  limit  determined  by  solution  freezing.  The 
lowest  temperatures  reached  were  210  K  for  the  30  and  40Z  NH^CHO  solutions.  The 
careful  sectional  baseline  deterainations  which  made  1Z  accuracy  aeasureaents 
possible  in  our  previous  work,^  were  again  utilised,  and  results  of  better 
than  1Z  precision  were  consistently  obtained. 

Expansivities  were  measured  in  the  stable  solution  temperature  region  by 
using  a  dilatoaeter  aade  of  precision  bore  capillary  (=  2  an  in  diaaeter)  with 
a  glass  bulb  attached  to  one  end.  The  conversion  of  the  meniscus  position, 
read  with  a  cathetoneter,  into  the  voluae  was  aade,  after  calibration  at  20°C 
with  carefully  degassed  aercury.  The  saall  corrections  for  pyrex  sxpansivity 
wss  aads  using  the  literature  (CXC  Handbook  of  Chemistry  sad  Physics)  value  of 
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2.5  x  10'  K~*.  The  expansivity  values  calculated  from  the  smoothed  volume  va. 

temperature  curve  were  obtained  with  an  uncertainty  of  ca.  +  1Z  reproducibly 
for  different  series  of  sample  preparation  and  measurement. 

Compressibility  measurements  were  begun  by  determining  absolute  values  for 
each  solution  at  5  MPa  and  25°C  using  a  30  cm  long  12  mm  precision  bore  pyrex 
capillary  sample  container  in  which  the  sample  was  isolated  from  the  pressure 
transmitting  fluid  by  a  5  cm  Hg  slug.  The  effect  of  pressure  on  the  internal 
volume  of  the  capillary  vessel  was  calibrated  by  making  similar  measurements  in 
the  same  vessel  on  spectroscopic  grade  benzene  to  obtain  the  glass  dilativity, 
d  [d  •  Cl/V  3 V/3F£tlterTial >  •  -  (compressibility)].  The  1  atm  pressure  volume 
of  the  vessel  was  determined  using  degassed  mercury.  To  make  high  pressure  and 

low  temperature  measurements  on  these  same  solutions,  the  fine  capillary 

(8) 

mercury  bead  displacement  method  of  Kanno  was  adapted  using  the 
above-measured  absolute  values  to  determine  tbe  effective  capillary  length. 
Rather  than  correct  results  for  the  capillary  dilation  effect  as  in  the  earlier 
measurements,  we  included  the  effect  of  dilation  in  the  length  calibration 
procedure  itself  by  least  squares  fitting  measured  bead  displacements  in  the 
range  0-10  MPa  to  the  expression 


3(i-£  ) 

l  d  +  (l-H  )d+ - 

o  o _ 3P 

kt  =  h  +  a-*.  > 

o  o 


(2) 


where  is  the  accepted  compressibility^*^  in  the  known  temperature  range, 
(£-£0)  is  the  difference  between  the  reading  of  the  reference  and  meniscus 
points,  and  d  is  the  glass  dilativity.  The  dilativity  thus  determined 
simultaneously  with  the  effective  length,  was  found  to  be  the  same,  3*67  x  10“® 
MPa-1  within  ♦  10Z  aa  that  obtained  in  the  bensene  calibration  series  with  the 
precision  bor*  capillary  of  diameter  1.0  ta.  An  accuracy  of  10Z  is  sufficient 
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to  reduce  the  error  in  compressibility  due  to  dilation  to  less  than  0.1Z. 

Values  of  <T  obtained  for  pure  water  in  the  range  0  to  -20°C  were  found  to  be 

(8) 

the  same  as  the  corrected  values  reported  earlier  to  within  the  IX 
uncertainties  of  the  measurement. 

Results 

Results  of  heat  capacity  measurements  on  the  binary  solutions  in  the 

temperature  interval  between  290  K  and  210  K  are  presented  in  Fig.  1. 

Measurements  close  to  pure  water  have  not  been  attempted  as  (a)  their  general 

( 2) 

form  has  been  shown  in  an  earlier  paper,  and  (b)  measurements  in  the 
supercooled  region,  requiring  emulsion  techniques,  are  rather  tedious. 

Comparison  of  our  solution  data  can  be  made  with  previous  work  by  Visser  et 
al*6^  at  298.15  K.  Over  most  of  the  composition  range  we  find  agreement  within 
IX,  which  is  the  uncertainty  of  our  measurements.  To  estimate  the  "normal" 
part  of  the  pure  water  heat  capacity  a  weakly  curvilinear  function  seems 
required  by  the  solution  data.  The  curvature  is  such  as  to  emphasise  the 
extraordinary  upswing  in  Cp  as  pure  water  is  approached  at  low  temperatures.  K 
crude  estimate  of  the  partial  molar  heat  capacity  of  RH^CHO  at  high  dilution  at 
-35°C  yields  the  extraordinary  value  -900  J  K-1  mol*"1,  comparable  with  the 
value  -850  J  K  *  mol  *  for  HjOj  quantitatively  obtained,  but  not  reported,  from 
our  earlier  emulsion  data  on  the  Hj02  ♦  HjO  system.* 

Results  of  the  expansivity  measurements  on  the  binary  solutions  are  shown 
in  Figs.  2  and  3.  Figure  2  shows  the  raw  data  for  each  composition  vs. 
temperature  to  indicate  the  precision  of  the  measurement.  Rote  that 
measurement  of  a  to  ♦  IX  requires  precision  in  volume  measurement  to  1  in  10^. 
Data  in  Fig.  3  are  obtained  from  best  lines  through  the  Fig.  2  points.  In  this 
case  the  relative  magnitude  of  the  low  temperature  anomaly  encountered  on 
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approaching  pure  water  is  very  large,  since  not  only  is  there  a  reversal  in 
sign,  but  the  absolute  magnitude  of  the  property  changes  fourfold  between  50 
and  100X  H^O.  As  in  the  case  of  the  solution  heat  capacity  data,  the  most 
plausible  extrapolation  of  the  binary  solution  range  is  a  curvilinear  one,  our 
best  representation  of  which  is  shown  by  the  dashed  lines  at  300,  270,  and 
240  K  in  Fig.  3.  In  the  case  of  the  expansivities,  there  are  no  literature 
data  with  which  comparisons  can  be  made. 

Results  of  the  compressibility  <T  measurements,  in  this  case  reported  for 
10  MPa  are  shewn  in  Fig.  3.  As  in  the  case  of  the  expansivity  data,  an 
influence  of  the  water-water  interactions  becomes  evident  at  much  larger  mole 
fractions  of  the  second  component  than  in  the  case  of  the  heat  capacity  data,  a 
point  which  is  discussed  farther  below.  Particularly  at  the  higher 
temperatures,  the  binary  solution  behavior  is  very  simple,  and  the 
extrapolations  to  0Z  formamide  content  are  fairly  unambiguous.  Again,  there 
are  no  literature  data  with  which  comparison  can  be  made. 

Values  of  are  also  available  for  high  pressures,  up  to  1500  bar  but  they 
will  be  considered  in  a  separate  paper. 

Discussion 

In  Figs.  5,  6  and  7  we  show  the  experimental  value  of  the  three  derivative 
properties  and  the  background  component  obtained  from  the  extrapolated  value  of 
the  property  at  CHQ  for  each  isotherm  of  Figs.  1-3. 

In  Fig.  5  we  include  the  background  Cp  functions  obtained  from  similar 

(23 

studies  on  the  HgOg  -  HgO,  and  ♦  HjO  systems.  The  present  system 
yields  values  equal  to  those  from  -  R^O  solutions  at  high  temperatures  and 

values  intermediate  between  and  HgOg-HgO-based  values  at  low 

temperatures  due  to  a  smaller  temperature  dependence  than  in  the  _  g^O 
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case.  This  is  also  Che  behavior  of  the  "normal"  component  extracted  from 
LiCl-HjO  solution  studies  to  be  reported  separately,  indicating  that  some 
average  of  the  values  observed  to  date  lying  closer  to  than  to 

based  values  is  the  appropriate  background. 

In  the  case  of  the  expansivity,  Fig.  6,  it  is  possible,  with  clarity,  to 

show  the  anomalous  part  of  the  expansivity  as  veil  as  the  total  and  background 

components  from  which  it  is  obtained  by  difference.  To  check  the  sensitivity 

of  our  treatment  to  the  element  of  arbitrariness  which  enters  in  our 

extrapolation  function,  we  show  "normal"  and  derived  "anomalous"  components  for 

three  different  choices  of  extrapolation  function.  This  is  a  more  serious 

problem  for  the  volume  derivative  properties  than  for  C^,  at  least  in  the 

present  system,  because  the  effect  of  the  water  structure  is  felt  out  to  ssich 

larger  mole  fractions  of  the  2nd  component  (compare  Figa.  1  and  2).  The  three 

backgrounds  in  Fig.  6  differ  in  the  number  of  mole  fraction  points  used  in  the 

quadratic  functions  which  were  extrapolated  to  X  ■  0.  Extrapolation  (a)  shown 

in  Fig.  3  used  1.0,  0.8,  and  0.6  mole  fraction  of  N^CHO,  extrapolation  (b) 

used  these  points  and  also  KRHjCHO)  ■  0.4,  while  extrapolation  (c)  used  a 

straight  line  thorugh  the  first  two  points  at  0.8  and  1.0  mole  fractions  of 

HH2CH0.  He  will  see  below  that  while  this  small  difference  has  little  effect 

on  the  best  fit  value  of  Tg  in  Eq.  (1)  it  has  a  significant  effect  on  the 

(3) 

associated  exponent  y.  The  background  suggested  by  Speedy  in  the  absence  of 
any  experimental  guidance,  is  identical  to  (c)  shown  as  a  dashed  line. 

In  Fig.  7,  we  show  the  total  and  background  components  of  the 
compressibility  obtained  by  the  straightforward  extrapolations  of  Fig.  4. 
Included  for  comparison  is  the  measured  compressibility  of  the  simplest  diol, 
ethylene  glycol,  one  of  the  few  other  simple  molecular  liquids  with  a 
compress iblity  smaller  than  that  of  water. 


A  problem  for  our  subsequent 
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analysis  is  seen  immediately  in  the  tendency  of  the  background  component  to 
become  equal  to  the  total  compressibility  i.e.  anomalies  vanish,  at 
temperatures  as  low  as  140  K . 

The  background  obtained  in  the  present  work  is  much  larger  than  that 

assumed  (by  relation  to  ethylene  glycol)  in  an  earlier  analysis  of  the 

(8) 

compressibility  anomaly,  of  water.  It  is  also  larger  than  the  "normal" 
component  for  water  recently  deduced  by  Conde  et  al^  from  sound  velocity 
studies  on  ethanol  +  water  solutions.  In  the  latter  system,  the  anomalous 
temperature  dependence  of  compressibility  is  eliminated  by  addition  of  just  10 
mol  X  KtOH,  whereas  obvious  anomalies  persist  even  at  35Z  formamide  in  the 
present  system.  The  molecular  shapes  and  siaea  of  the  two  second  components 
are  rather  similar  (65^ — and  the 

difference  in  behavior  must  reflect  the  influence  of  the  hydrogen  bonding  amine 
group  in  formamide  in  increasing  the  structural  compatibility  with  water.  The 
suggestion  is  that  formamide  may  be  too  compatible  with  water  for  our  purpose 
of  separating  out  abnormal  fluctuation  components,  and  that  the  background 
extracted  in  our  work  is  therefore  an  overestimate.  The  contrast  serves  to 
illustrate  the  severe  limitations  of  the  binary  solution  extrapolation  method 
for  separating  "anomalous"  from  "background"  contributions  to  an  intensive 
thermodynamic  property.  One  runs  the  danger  of  choosing,  as  an  "appropriate” 
binary  system  for  background  determination  the  one  which  finally  delivers  the 
"anomalous”  component  one  wants  to  find. 

With  this  chastening  thought  in  mind,  let  us  observe  the  temperature 
dependence  of  the  "anomalous"  components  of  water  for  the  three  properties  Cp, 
a  and 

We  show  graphically  the  tests  of  Eq.  (1)  for  and  a,  respectively,  in 


Tigs.  8  and  9 
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The  corresponding  plots  for  the  excess  <T,  which  cover  the  smallest 
supercooling  range  (T  >  -26°C),  are  not  shown  because  the  parameters  which  give 
the  best  linear  fits  are  absurd  -  for  instance,  Tg  -  -100  K.  Although  this 
sort  of  result  would  seem  to  immediately  invalidate  our  approach  we  must  note 
that  the  excess  obtained  by  Conde  et  al^  from  a  similar  analysis  of  data 
obtained  on  ethanol  ♦  H^O  solutions  (from  light-scattering  measurements)  fitted 
Eq.  (1)  with  parameters  very  similar  to  those  of  our  present  best  fits  to  the 
excess  and  a  values,  vis.^^  Tg  “  228  K,  Y  ■  0.91.  This  observation 
encourages  us  to  look  more  closely  at  the  dependence  of  the  Iq.  (1)  parameters 
for  the  excess  a  data  on  the  composition  range  of  binary  data  used  to  determine 
the  background  component.  The  point  is  that  C^,  which  is  determined  by  entropy 
fluctuations  (tt(AS)  )  seems  to  become  "normal"  in  solutions  of  quite  small  mole 
fractions  of  added  formamide  (Pig.  2),  while  k^,  which  is  determined  by  volume 
fluctuations  alone  («(AV)2),  may  be  "abnormal"  across  the  whole  composition 
range,  a,  which  is  determined  by  the  product  of  density  and  volume 
fluctuations,  might  therefore  be  expected  to  show  some  diagnostically  useful 
trend  as  the  composition  range  from  which  the  background  contribution  is 
assessed,  is  changed.  In  Fig.  7  the  three  best  fit  lines  (a),  (b)  and  (c),  are 
for  the  excess  (anosMlous)  components  of  a  obtained  when  successively  smaller 
composition  ranges  (which  successively  weaken  any  residual  anomalous  water 
influence)  are  used  in  the  background  determination  (see  previous  section). 

Apart  from  further  demonstrating  the  weakness  of  our  binary  solution  method 
for  separation  of  "normal"  and  "anomalous"  components  of  water  properties,  this 
sequence  is  interesting  in  that  it  shows  the  best  fit  values  of  y^  decreasing 
from  high  values  to  unity  as  the  possibility  of  residual  cooperative  water 
influence  is  increasingly  reduced. 

Although  this  might  seem  to  reinforce  the  significance  of  the  value  1.0  for  y 
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suggested  in  our  earlier  analyses ,  an  extension  of  the  line  of  reasoning  would 
imply  s  smaller  value  of  Y  as  a  possibility  to  be  dealt  with.  If  the 
persistence  of  abnormal  volume  fluctuations  across  the  binary  is  so  influential 
then  Cp  data  should  yield  the  most  reliable  separations,  hence  the  most 
reliable  values  of  Y.  Cp  data  are  also  available  over  the  widest  temperature 
range  and,  perhaps  marginally,  with  the  greatest  accuracy  (we  have  not  used 
data  below  -35°C  because  of  the  difficulty  of  proving  absence  of  incipient 
crystallization  distortions  of  the  data).  The  present  Cp  data  yield  "normal" 
component  magnitudes  very  similar  to  those  obtained  from  an  earlier  ♦  HjO 

solution  study,  but  with  a  smaller  temperature  dependence  (see  Fig.  5).  These 
distinctions  lead  to  a  value  of  Y  ■  0.75  which  is  the  smallest  yet  obtained. 
However,  such  a  small  value  of  Y  forces  T^  to  the  correspondingly  large  value 
of  210  K  which  seems  improbable.  For  the  H^H^-based  Cp  background,  we  earlier 
obtained  Y  “  0.91,  Ta  -  226  K.  Without  knowledge  of  the  volume  derivative 
properties  based  on  as  a  second  component  we  cannot  comnt  on  which  if 


either  set  to  favor. 

It  seems  reasonable  overall  to  conclude  that  if  Eq.  (1)  is  to  apply  to 
these  phenomena  then  Y  values  close  to,  though  less  than,  1.8  seem  the  most 
plausible. 

Values  close  to  unity  would  be  in  accord  with  our  earlier  understanding  of 

mean  field  theory  predictions  for  approach  to  a  mechanical  instability  for  a 

f  1 2) 

liquid  (e.g.  on  superheating),  '  but  it  now  appears,  from  analysis  by 

Speedy ^ ^  and  earlier  work  by  Compagner/*^  that  series  expansions  around  a 

spinodal  point  actually  require  Yc  -  Ya  -  yk  -  0.5.  The  exponent  1.0  only 

p  T 

applies  at  the  critical  density  and  leads  to  thermodynamic  inconsistencies  at 
T#  for  other  densities.  It  is  not  clear,  however,  how  close  to  T#  the 
observations  must  be  made  for  the  exponent  0.5  to  be  observed.  Ve  recently 
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noted^  that  the  Haar-Gallagher  equation  of  atate  for  water^^  gives  behavior 
near  the  superheated  liquid  spinodal  point  which  is  compatible  with  an  exponent 
of  unity  for  values  in  the  range  of  data  for  supercooled  water,  but  went 
over  to  Y  ~  0.5  behavior  when  the  actual  instability  temperature,  and  data  very 
close  to  it,  were  used  in  the  analysis.  So  values  of  y  near  unity  from  the 
present  study  are  not  incompatible  with  a  spinodal  interpretation  of  the 
observed  anomalies.  They  are  also  compatible  with  the  alternative 
possibility^7^  that  the  anomalies  reflect  the  approach  to  a  higher  order 
transition  within  the  liquid  state  associated  with  the  establishment  of  some 
long-range  order  in  the  hydrogen  bond  network  -  a  notion  which  has  recently 
been  provided  with  a  topological  basis  by  Stanley  and  Teixeira.^ 

Thus  the  issue  of  the  origin  of  the  curious  behavior  of  water  at  low 
temperatures  cannot  be  settled  by  these  measurements.  The  difference  between 
the  apparently  distinct  postulates  of  spinodal  (edge  of  liquid  free  energy 
surface)  and  internal  transition  (bump  on  the  liquid  free  energy  surface)  may 
not  be  as  great  as  it  seems.  Indeed  it  is  not  inconceivable  that  a  changeover 
could  occur  from  one  form  of  behavior  to  another  as  pressure  is  increased 
sufficiently  above  1  atm.  One  published  equation  of  state^^  hints  at  this 

•  •  •  ( to) 

possibility  in  its  predictions  of  properties  well  outside  its  valid  range. 

Such  oddities  would  be  a  property  of  the  liquid  free  energy  surface  and  would 
have  nothing  to  do  with  the  equilibrium  between  water  and  ice  which  can  be 
observed  in  separate  experiments. 

Anomalies  at  Constant  Volume 

Since  the  extraordinary  nature  of  low  temperature  water  is  clearly 
dependent  on  the  fluctuations  in  entropy  and  volume  made  possible  by  the 
migration  of  hydrogen  bonds  between  water  molecules  which  continuously  attempt 
to  establish  a  local  open  tetrahedral  network  (amorphous  ice-like)  order, 


certain  difference*  in  the  temperature -dependent  properties  from  thoae 

diacuaaed  above  muat  be  expected  when  measurements  are  performed  at  eonatant 

volume  or  constant  entropy.  Indeed,  the  vanishing  of  the  heat  capacity 

divergence  and  its  replacement  with  a  weak  maximum  when  examined  at  constant 

volume  was  reported  at  the  beginning  of  the  current  wave  of  supercooled  water 
(19) 

Studies.  Likewise,  evidence  for  a  maximum  in<f  has  recently  been  seen  in 

the  clever  levitated  droplet  sound  velocity  measurements  of  Trinh  and 
(20) 

Apfel.  The  latter  maximum  has  also  been  implied  by  short  extrapolation  of 

similar  data  from  recent  light  scattering  studies  on  EtOH  +  HjO  solutions  in 
which  the  k  maxima  were  directly  observed/7^  It  is  therefore  of  interest  to 
use  the  present  data  to  see  whether  there  is  an  important  "anomalous  component" 
to  the  heat  capacity  for  water  when  studied  at  constant  volume. 

in  the  range  17  to  -15°C  can  be  obtained  from  the  present  data  using  the 
standard  thermodynamic  relation 

Cy  -  Cp  -  VT  «2/<t  (3) 

The  results  are  shown  in  Fig.  10.  Before  commenting  on  their  implications  it 

should  be  remembered  that  the  values  of  C  thus  calculated  are  (3E/3T)  at  the 

v  v 

volume  jyt  which  is.  measured .  rather  than  at  a  single  value  of  volume  which 
is  the  same  for  all  temperatures.  For  our  main  purpose  of  observing  the 
extrapolation  of  solution  values  to  pure  water,  however,  this  distinction  is 
iaaaterial,  though  it  affects  the  form  of  vs.  T  shown  in  the  inset  to  Fig. 
10.  At  a  single  fixed  volume  the  vs  T  function  will  be  somewhat  sharper, 
depending  on  the  value  of  volume  held  constant. 

Figure  10  clearly  shows  that  there  is  indeed  a  large  contribution  (*  12X) 
to  0^  of  pure  water  from  an  aspect  of  water-water  interactions  which  is 


n 

eliminated  by  addition  of  only  20  mol  Z  of  formamide. 

This  ia  not  surprising:  the  density  and  temperature  fluctuations  will 

occur  whether  or  not  the  system  can  change  it  volume.  It  is  useful,  however 

because  it  allows  us  to  form  a  clearer  idea  of  the  magnitude  of  the 

specifically  "waterlike"  fluctuations  which  determine  for  water,  and  to 

obtain  a  form  for  the  "excess  C  "  as  a  function  of  T  which  is  useful  for 

v 

theoretical  purposes  (most  theories  being  constructed  for  constant  volume 
systems).  and  for  pure  water,  are  shown  against  the  "normal" 
component  obtianed  by  the  solution  extrapolations  in  the  upper  insert  to  Fig. 
10. 

The  values  of  shown  in  Fig.  10  are  assessed  using  the  most  recent 

measurements  of  for  supercooled  water^  which  revise  downward  by  some  5Z  at 

( 19) 

the  lowest  temperature,  the  values  first  reported  for  this  property.  The 

revision  has  the  effect  of  producing  a  maximum  and  a  sharp  fall-off  in  at 

lower  temperatures,  where  previously  only  a  weak  maximum  had  been  seen.  The 
point  at  -35°C  (238  K)  is  quite  uncertain,  since  it  depends  on  combinations  of 
measurements  of  intensive  thermodynamic  properties,  each  of  which  is  difficult 
to  measure  with  accuracy,  and  all  of  which  are  changing  very  rapidly  with 
temperature  in  this  region.  Furthermore,  the  value  of  used  for  this  point 
required  a  9°C  extrapolation  of  existing  data.  It  seems  very  improbable  tbat 
Cv  for  pure  water  would  fall  below  the  value  of  the  "normal"  component  assessed 
in  the  present  work,  hence  we  suggest  that  the  most  probable  form  for  is 
tbat  shown  by  the  dashed  line,  which  approaches  the  (normal)  line 
asymptotically. 

The  separated  excess  due  to  that  component  of  the  mean  square 

temperature  fluctuation  that  is  uniquely  "water-like”,  is  shown  as  AC  in  the 

v 

lower  insrt  to  Fig.  10.  It  has  the  form  of  a  Shottky  anomaly  which  is 
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predicted  by  the  simple  random  network  "bond  lattice"  model,  in  which 

H -bonds  forming  during  thermal  de-excitation  are  assumed,  in  the 

non-cooperative  case,  to  cause  the  same  volume  and  energy  increments 

irrespective  of  their  number.  In  this  it  differs  from  the  Stanl ey-Teixeira 

random  bond  percolation  model ^  which  aasigna,  very  reasonably,  an  extra 

volume  per  bond  to  the  "saturated"  cluaters  of  bonds  which  arise  naturally  in 

the  course  of  cooling  at  constant  pressure.  A  similar  idea  waa  briefly 

described  in  an  earlier  paper  in  which  it  waa  pointed  out  that  a 

(  221 

volume-cooperative  modification  of  the  bond  lattice  model  (in  which  the  AV 

per  bond  is  made  proportional  to  the  mole  fraction  of  "on”  bonds)  could 

qualitatively  describe  the  anomalous  compressibility  observations  for  water. 

Since  according  to  Le  Chatelier's  principle  the  formation  of  large  volume  bond 

clusters  ( "patches") ^  natural  to  random  bonding  systema  of  HjO  molecules^ 

at  low  temperatures  must  be  discouraged  by  a  constant  volume  condition,  fixing 

the  system  volume  must  introduce  a  non-random  element  into  the  bond 

distribution.  This  is,  of  course,  manifested  as  an  increase  in  pressure  with 

decreasing  temperature  at  constant  volume. 

The  form  of  the  anomaly  is  not  unlike  that  of  the  adiabatic 

compressibility  implied  by  the  sound  velocity  studies  of  Trinh  and  Apfel,^^ 

and  of  Conde  et  al^  referred  to  above,  and  the  temperatures  of  their  maxima 

differ  by  only  some  10°.  in  the  RgO-etbanol  solutions  of  the  latter  study  the 

<s  maximum  is  sharply  defined.  Temperature  and  volume  fluctuations,  which 

(21) 

determine  and  <a,  are  statistically  independent ,  however,  and  it  is 

(1)  2 
worth  remember ing  that  the  isothermal  compressibility  (*  (A?)  )  increases 

4 

almost  as  rapidly  with  decreasing  temperature  at  a  constant  volume  of  18  cm 
mol-1  aa  it  does  at  a  constant  pressure  of  1  atmosphere. 
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Figure  Captions 

Figure  1  Expansion  coefficients  for  water-foraaaide  solutions  as  a  function 
of  temperature.  Each  point  is  deterained  from  coluan  length 
measureaents  taken  =  2.5  K  apart  in  teaperature.  Points  at 


*iih2cho  “  0  .40  represent  two  different  runs, 
scatter  at  this  composition  than  at  others. 


there  being  aore 
(Note  +  IX  indicator 


bars.) 


Figure  2  Heat  capacities  of  water-foraamide  solutions  at  temperatures 

between  210  and  290  K,  obtained  froa  smoothed  plots  of  raw  C  vs. 

P 

T  data  at  each  composition  indicated.  Values  may  be  compared  with 
those  of  Visser  et  al^)  at  298.15  K.  Connections  between  the  20 
mol  X  solution  and  pure  water  data  are  indicated  by  dashed  lines 
(for  measurements  in  this  difficult  region  see  ref.  2).  Natural 
extrapolations  of  solution  data  to  OX  foraaaide  are  shown  as  full 
lines.  Intercepts  at  X  ■  0  indicate  "normal "  component  of  Cp  for 
pure  water. 

Figure  3  Expansion  coefficients  for  wster-formamide  solutions  at 

temperatures  between  240  and  300  K.  Alternative  possible 
extrapolations  to  obtain  the  "normal"  component  of  the  water 
expansivity  are  indicated  by  dashed  lines  (see  text  for  discussion 
of  these  alternatives). 

Figure  4  Isothermal  compressibilities,  at  the  pressure  of  10  KPa  for 

vater-foraanide  solutions  at  temperatures  between  233  K  and  293  K. 
Dashed  lines  indicate  plausible  extrapolations  to  obtain  a 
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Figure  5 


Figure  6 


Figure  7 


Figure  8 


"norMl"  component  of  the  compressibility  of  pure  water. 

The  "normal"  coaponent  of  the  heat  capacity,  Cp  of  water  obtained 

from  extrapolations  of  Fig.  2  for  the  present  water-foraamide 

solution  study.  Comparisons  are  made  with  "normal"  components 

( 2) 

obtained  previously  from  HjOj-H.O  and  data. 

The  figure  compares  these  background  components  with  the  total 
for  water  (see  ref.  7). 

Plots  of  the  "normal",  "anomalous"  and  total  (experimental) 
expanaivities  of  water  as  a  function  of  temperature.  Curves  (a), 
(b)  and  (c)  indicate  the  consequences  of  different  choices  of 
binary  solution  composition  range  for  extrapolations  to  0  water 
content  (see  Fig.  3).  "Anomalous"  component  suggested  in  ref.  3, 
is  shown  as  a  dashed  line. 

Comparison  of  the  total  isothermal  compressibility  of  water  with 
the  "normal"  component  of  k ^  extracted  from  the  present 
water-foraamide  solution  data,  and  from  the  ethanol-water  solution 
data  (of  ref.  1). 

Best  fit  of  the  "anomalous  component"  of  C  for  water  to  Eq.  (1) 

P 

(curve  a)  and  tests  of  Bq.  (1)  descriptions  of  the  "anomalous” 
component  of  C^  for  pure  water.  Curve  (a)  is  an  overall  best  fit, 
while  curve  (b)  is  a  constrained  best  fit  for  an  assigned  value  of 

y*  y  -  l. 
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Figure  9  (a)  Overall  best  fits  of  the  "anomalous"  components  of  the 

expansion  coefficient  of  pure  water,  using  the  three  alternative 
"background"  or  "normal"  component  extrapolations  discussed  in  the 
text.  Comparisons  indicate  sensitivity  of  the  y  parameter  to  the 
"background"  form. 

(b)  Constrained,  ym  1,  best  fits  of  anomalous  component  of  the 
expansivity  for  the  three  different  choices  of  background. 


Figure  10 


Plot  of  the  constant  volume  heat  capacity  for  vater-formamide 
solutions  assessed  from  combined  data  of  this  study. 


Extrapolation,  to 
for  water. 


0,  yield  a  "normal"  component  of  Cy 


Insert  1 :  "normal"  component  of  Cy  for  pure  water  compared  with 
total  Cy  and  total  Cp  for  water.  The  lowest  temperature  Cv  point 
is  very  uncertain,  and  low  temperature  form  of  C^  is  probably 
asymptotic  to  the  "normal"  component.  The  croases  are  Cy  values 
calculated  for  a  constant  value  of  V  -  18.0  cc/mole. 

Insert  2;  "Anomalous"  component  of  Cy  of  water  obtained  from 
insert  1  by  subtraction.  Note  Shottky  anomaly  form  for  this 
component . 
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